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Abstract 

A search has been performed for interstellar CH3D via its Jk = lo - Oo 

transition at 230 GHz and its Jk = 2o - 1 0 and Jk = 2i - 1 1 lines at 465 GHz using 

the NR AO 12 m and CSO 10 m telescopes towards Orion-KL. This search was 

done in conjunction with laboratory measurements of all three transitions of CH3D 

using mm/sub-mm direct absorption spectroscopy. The molecule was not detected 

* 

down to a 3a level of Ta < 0.05 K towards Orion, which suggests an upper limit 
to the CH3D column density of N < 6 x 10 18 cm- 2 in the hot core region and a 
fractional abundance (with respect to H2) of < 6 x 10"6 These measurements 
suggest that the methane abundance in the Orion hot core is f < 6 x 10’ 4 , assuming 
D/H ~ 0.01. Such findings are in agreement with recent hot core chemical models, 
which suggest CH4/H2 ~ 10 -4 . 

Subject Headings: ISM: abundances, ISM: clouds, ISM: molecules - molecular 
processes 
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1. Introduction 

Chemical models predict that CH 4 , the simplest stable organic molecule, to 
be one of the most abundant polyatomic species in interstellar dense clouds, since 
it only involves one carbon atom bonded to four hydrogen atoms (e.g. Mitchell 
1977, Brown & Rice 1986). In addition, methane is thought to have an enhanced 
abundance in hot, dense gas, such as the Orion hot core region due to evaporation 
of this molecule from grain mantles (e.g. Brown, Chamley, & Millar 1988). 
However, CH 4 is completely symmetric, having the geometry of a perfect 
tetrahedron, and hence does not possess a permanent electric dipole moment. 

Thus, the molecule does not have a normal pure rotational spectrum and 
consequently cannot be readily observed in cooler, dense gas typical of molecular 
clouds. 

Methane, on the other hand, does possess a rovibrational spectrum in the 
infrared which has been detected towards IRC+10216 (Hall and Ridgway 1978) 
and possibly in a few star-forming regions (Lacy et al. 1991). However, the kinetic 
temperatures found in molecular clouds, even toward star-forming regions, are 
typically only T ~ 10-100 K. Therefore, vibrational excitation is not commonly 
found in interstellar molecules, especially for regions of extended gas, the 
dominant constituent of molecular clouds. Hence, it is difficult to test models of 
interstellar chemistry using CH 4 infrared observations. 

An alternative method to obtain interstellar methane abundances is possibly 
by observing its singly-deuterated form, CH 3 D, which has a pure rotational 
spectrum. Because CH 3 D is only slightly asymmetric, however, its dipole moment 
is small: ji ~ 0.006 D (Wofsy, Muenter, & Klemperer 1970). Moreover, given the 
normal cosmic D/H ratio of 10‘ 5 , detection of CH 3 D should prove difficult. A 
search was carried out for CH 3 D via its Jk = lo - Oo transition at 232 GHz by 
Pickett, Cohen, and Phillips (1980), but given the receiver technology at the time, 
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not very sensitive limits (Ta* ~ 0.2 K) were obtained in Orion-KL. However, 
recent observations of interstellar molecules have clearly demonstrated that 
significant deuterium fractionation may occur in dense clouds and that the cosmic 
D/H ratio has little to do with abundances of deuterated species in comparison with 
their parent molecules. For example, considerable enrichment is found in several 
deuterated molecules toward Orion-KL, with HDCO/H 2 CO ~ 0.01-0.03, 
CH 3 OD/CH 3 OH ~ 0.01-0.06, and NH 2 D/NH 3 ~ 0.003-0.05 (Loren and Wootten 

1985; Mauersberger et al. 1988; Walmsley et al. 1987). 

Here we present results of a renewed search for interstellar CH 3 D towards 
Orion-KL. We have conducted measurements of the Jk = lo - Oo transition at 232 
GHz and the Jk = 2i - 1 1 and 2o - lo lines at 465 GHz of this isotopomer. We also 
have carried out laboratory measurements of all three CH 3 D transitions. Although 
spectral lines were detected at the Jk = 1 — » 0 and Jk = 2o - lo frequencies in 
Orion, no features were present at the Jk = 2i - 1 1 line, hence making our search 
negative. In this paper both the astronomical and laboratory data are summarized. 
We also present upper limits for the methane abundance and discuss its 
implications for interstellar chemistry. 

2 . Observations 

The fundamental transition Jk = lo - Oo of CH 3 D at 232,463.3 MHz (see 
Pickett, Cohen, and Phillips 1980) was originally searched for using the NRAO 1 
12 m telescope at Kitt Peak, AZ in April 1991. At this frequency, the beam size 
was 28” and the beam efficiency was T|b = 0-4. The temperature scale was 
determined by the chopper wheel method, corrected for forward spillover losses, 


1 NRAO is operated by the Associated Universities, Inc., under contract with the National 
Science Foundation. 
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and is given in terms of Tr*. Conversion to radiation temperature is then Tr = 
Tr*/tjb- The receiver used was a dual channel SIS mixer. The backends used 
were two 256 channel filter banks with 1 MHz resolution, one for each receiver 
channel. Data was taken in a double sideband (DSB) mode. 

Additional observations were done of the Jk = lo - Oo line using the Caltech 
Submm Observatory (CSO) in 1993 September. At 232 GHz, 85 ~ 40" and T|b = 

0.8. The receiver used was a single channel SIS mixer, again operated in DSB 

mode. In addition, measurements of the two K components of CH3D near 465 

GHz were conducted 1992 March and 1995 February, again using the CSO. The 

two K-components were observed simultaneously. The beam size at this frequency 

is 0b ~ 15" and t|b = 0.5. Again, the receiver used was a single channel SIS mixer. 

The temperature scale for all CSO observations is in terms of Ta*, such that Tr = 

* 

Ta /t|b. The backend used was a 1024 channel AOS spectrometer with 500 kHz 

resolution. All observations were carried out towards the position a = 5 h 32 m 46 s 8; 
5 = -5°24’23".0 (1950.0) in Orion. 

3. Experimental 

Spectra of CH3D were obtained in the laboratory using a millimeter/sub- 
milliliter wavelength direct absorption spectrometer, which is described in detail 
elsewhere (Ziurys et al. 1994). Briefly, the instrument consists of a tunable source 
of mm/sub-mm radiation, a gas absorption cell, and a detector. The source for the 
spectrometer are phase-locked Gunn oscillators (65-140 GHz), used in conjunction 
with frequency multipliers to cover the range 1 15-520 GHz. The radiation is 
quasi-optically propagated through the absoiption cell, which is a double-pass 
system. The detector is a helium-cooled InSb bolometer. 

For these experiments, CH3D gas, purchased from a commercial source 

(Icon), was sealed in a quartz glass tube about 6 cm in diameter and 400 m in 
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length with quartz flats sealing both tube ends. CH3D was added to the cell at the 
desired pressure (~50 mtorr) and the cell then sealed. The tube was then inserted 
into the spectrometer and data taken, allowing for long signal averaging without 
having to free-flow CH3D gas. 

4. Results 

In the initial search for CH3D using the NRAO 12 m telescope, a weak 
feature was found at the frequency of the Jk = lo - Oo line at 232 GHz towards 
Orion-KL. This line was subsequently confirmed using the CSO, and this 
spectrum is shown in the top panel of Fig. 1. As can be calculated from Table 1, 
this feature has a linewidth of ~4 km/s and an LSR velocity near 5 km/s, if it arises 
from the Jk = lo - Oo line of CH3D. The frequency of this transition of CH3D had 
been previously measured in the laboratory by Pickett, Cohen, and Phillips (1980) 
to be 323,644.327(18) MHz, and hence there was no question as to its exact 
value. Such line parameters suggested that the observed feature, if due to CH3D, 
had an origin in the Orion hot core because of its slightly lower velocity (5 vs. 9 
km/s). Because large enhancements for deuterated species are commonly found in 
the hot core (e.g. Plambeck and Wright 1987), this result was not unexpected. A 
feature was also detected at near the CH3D frequency in SgrB2 and W51. 

The results of the search for the Jk = 2i - 1 1 and 2o - lo transitions near 465 
GHz were not as clear. In the March 1992 measurements, the data showed two 
spectral features which possibly could correspond to the Jk = 2i - 1 1 and Jk = 2o - 
lo components of CH3D, provided they had velocities near 5-6 km/s. However the 
candidate Jk = 2i - 1 1 feature appeared as a shoulder on a stronger line, which is 
unidentified, and the signal-to-noise was not particularly good. Moreover, the 
frequencies of these two K components of the J = 2 — » 1 transition of CH3D had 
never been directly measured before in the laboratory, and could only be predicted 
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from the ro-vibrational data of Tarrago et al. (1976). The absolute accuracy of 
such calculations was not clear. 

Because of these uncertainties, we decided to directly measure the J = 2 — » 1 
CH 3 D transitions in the laboratory, and remeasure the J = 1 — > 0 line frequency as 
well. The details of these measurements have been described in the experimental 
section of this paper. The resulting rest frequencies are listed in Table 2, and have 
an accuracy of better than ±75 kHz. A sample spectrum of the J = 2 — > 1 
transitions is shown in Fig. 2. 

With these new laboratory frequencies, we repeated our search for the 
J = 2 — > 1 lines of CH 3 D in the February 1995. For these observations, a spectrum 
with much better sensitivity was obtained towards Orion because of an improved 
receiver. These data are shown in the bottom panel of Fig. 1 . As this figure 
illustrates, a line with Ta* ~ 0.5 K is present at the frequency of the Jk = 2o - lo 
transition, for a velocity of Vlsr km/s, the same LSR velocity as the J = 1 — > 0 
candidate line. However, there is no obvious feature present at the frequency of 
the Jk = 2i - 1 1 transition, which is 15.15 MHz lower in frequency (or to the left) 
of the Jk = 2o - lo component. There are lines near 465,238 MHz and 465,232 
MHz, either which could correspond to the Jk = 2i - 1 1 line, but neither have a 
matching Jk = 2o - lo component ~15 MHz to higher frequency. 

The obvious failure to detect two distinct K-components of CH 3 D is clear 

evidence for its absence in Orion. The components are separated in energy by only 

8 K, so both should be detectable in the high excitation gas of the Orion hot core. 

An upper limit to the CH 3 D antenna temperature can be derived from the Jk = 2j - 

* 

1 1 component, which is Ta ~ 0.05 K. 

The frequencies, intensities, and linewidths of the spectral lines detected in 
Figure 1 are listed in Table 1. In addition, revised rotational constants for CH 3 D 
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are given in Table 2. As the table shows, they are in good agreement with past 
estimates from Pickett et al. (1980) and Tarrago et al. (1976). 

5. Discussion 

5. 1 Column Density and Fractional Abundance for CH3D 

The upper limit to the CH3D column density was calculated using the 
following formula, which assumes low optical depth, for (J + 1)k ~ * JK: 

3kl0 5 T R Av , /2 J 'e wtT “Q„ 
tol 87t 3 vp.J(J 2 - K 2 )e -AE/kTro, S 1K 

In this equation, which does not assume the Rayleigh-Jeans approximation, v is the 
frequency of the transition, |i 0 the dipole moment, AE the energy of the J^ 1 level 
above ground state, and Trot and T ex the rotational and excitation temperatures, 
respectively. Tr is the upper limit to the measured line temperature, and Avj /2 the 
assumed linewidth. The term Q rot is the rotational partition function, which for a 
symmetric top molecule can be approximated by (e.g. Townes and Schawlow 
1975): 


Q m ~Vb 2 C/*T 


( 2 ), 


where B and C are the rotational constants of the molecule. Sj k defines the 
statistical weight factor due to the presence of three equivalent nuclei with I = 1/2, 
i.e. the three protons (also see Townes and Schawlow 1975): 

2(4I 2 +41) 3 . 

4I 2 +41 + 1 
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In this formula, the 2.7 K microwave background was neglected (c.f. 

Ziurys, Hollis, and Snyder 1994). Also, because the dipole moment of CH3D 

is so small, it was assumed that T ex = T rot , i.e., the excitation temperature within 

the 2 i-l 1 transition equals that governing the population in the rotational ladder. 

To actually calculate the limit to the column density, it was assumed that 

CH3D would arise primarily from the Orion hot core region. Simple, heavily 

saturated species such as NH3 appear to high abundances in the hot core, and CH4 

is likely to follow this pattern (e.g. Brown, Chamley and Millar 1988). For this 

region, the approximate kinetic temperature is ~ 200 K (e.g. Blake et al. 1987), and 

hence it was assumed that T ex = T rot = 200 K. Also, Avi/2 was estimated to be 10 

km/s, the canonical “hot core” linewidth. The upper limit to the line radiation 

* * 

temperature was Tr= Ta /Tlg= 0-1 Km. where Ta < 0.05 was the upper limit for 
the Jk = 2i -1 1 transition. 

Using these assumptions, the upper limit to the CH3D column density in the 
Orion hot core was calculated to be N tot ~ 2 x 10 i8 cm -2 . This number is based 
solely on the limit obtained for the Jk = - 1 \ transition (assuming V^sr - 5 

km/s), which was not contaminated by other lines as the other two transitions were. 
This limit assumes the source fills the 15” beam at 465 GHz. Correcting for the 
hot core size of 10”, (Masson et al. 1985), however, increases the upper limit by a 
factor of 3, or N tot < 6 x 10 18 cm' 2 . If the total hydrogen column density is 
assumed to be N(H2) ~ 10 24 cm' 2 (e.g. Masson et al. 1985), as is appropriate for 
the hot core, then the fractional abundance of CH3D in this region is f ~ 6 x 10' 6 . 

5.2 Implications for the CH 4 Abundance. 

Estimating the methane concentration from CH3D depends critically on the 
D/H ratio, which fortunately, can be estimated from other molecules in Orion. 
Large deuterium enhancements have been found in several clumps towards the 
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KL/IRc2 region, including the hot core, the so-called “compact ridge”, and in the 
“northern condensation” (eg. Mangum, Plambeck, and Wootten 1991). In the hot 
core, measurements of NH2D/NH3 by Walmsley et al. (1987) yield D/H ~ 0.003, 

while DCN/HCN suggest D/H ~ 0.005 (Mangum, Plambeck, and Wootten 1991). 
Observations of CH3OH and CH3OD by Mauersberger et al. (1988) indicate higher 

ratios of D/H ~ 0.01-0.06. However, it is not entirely clear from what region the 
CH3OD lines actually arise. If a D/H ratio of 0.01 is assumed, the upper limit to 

the fractional abundance of methane in the Orion hot core is f ~ 6 x 10 -4 , with a 
corresponding column density limit of N tot (CH4) ~ 6 x lO 2 ^ cm' 2 . 

The only source where methane has conclusively been detected is in the 
circumstellar shell of IRC+10216 (Hall and Ridgway 1978). These authors 
observed the 3.3 pm V3 vibrational band of CH4, and derived a corresponding 
column density of N tot = 3 x IOI 7 cm* 2 . They also measured a CO column depth 
of Ntot = 10 20 cm' 2 in IRC+10216. If CO/H2 ~ 10* 4 , then the fractional abundance 
of methane, relative to H2, is f ~ 3x 10" 7 in this object. This value is several orders 
of magnitude below our upper limit of f ~ 6 x 10' 4 in Orion. However, the 
chemistry in the envelope of IRC+10216, a late-type carbon star, is quite different 
from that of the Orion hot core, so methane abundance are likely to vary between 
the two sources. 

Lacy et al. (1991) searched for the 7.6 |im V4 band of CH4 towards several 
molecular clouds. There authors appeared to detect the R(0) and R(2) lines of this 
band towards NGC7538 IRS9, although the spectra were highly contaminated by 
telluric features. From these two lines, Lacy et al. derived a CH4 column density 
of N tot ~ 2 x 1016 cm -2 towards NGC7538. Using their quoted CO column depth 
of 1.5 x 1019 cm -2 , and assuming CO/H2 ~ 10‘ 4 , the fractional abundance of CH4 
in this molecular cloud is f ~ 10' 7 , relative to H2. A similar fractional abundance 
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was possibly found for W33 A. Again, these CH4 abundances are several orders of 
magnitude lower than our upper limit. 

Lacy et al. also suggest a possible detection of the R(0) line of methane in 
Orion/IRc2. However, they do not derive a column density for CH4 from this 

measurement, although a comparison of methane abundances may have been 
useful. Lacy et al. determined an LSR velocity of the possible R(0) feature of 
Vlsr = -1 km/s. To our knowledge there is no known molecular emission in 

Orion that occurs at that velocity. 

5.3 Constraints for Chemical Models 

The “hot core” cloud model by Brown, Chamley, and Millar (1988) predict 
high abundances of simple saturated molecules such as NH3 and CH4. These 

calculations basically considered gas phase reactions and grain mantle evaporation 
in dense gas heated by a young star, with T ~ 150 K. In this model, these authors 
estimate f (CH4) ~ Kb 4 , relative to H2. A newer set of calculations by Brown and 

Millar (1989) basically consider the same model as Brown et al. (1988), but 
additionally consider some grain-surface reactions. Once again, the estimated 
abundance of methane for hot core-type gas is f ~ 10* 4 . 

The upper limit for methane in the Orion hot core region determined from 
the CH3D search suggests f ~ 6 x lO 4 . This value certainly does not contradict 

model predictions; however, it is not a stringent enough limit to put the chemical 
calculations to the test. Given the contaminating spectral features at two of the 
three easily accessible CH3D rotational transitions, obtaining more sensitive limits 

will be difficult. 

6 . Conclusions 

A search has been carried out for the = 1 q - 0q, 2q - lo» and 2j - 1 j, 
mm/sub mm transitions of CH3D at 232.64, 465.25, and 465.24 GHz. Laboratory 
measurement of these CH3D frequencies were also performed to substantiate 
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the astronomical work. Although spectral lines were clearly present at the lo - Oo 
and 2 q - lo frequencies for CH3D in Orion, assuming the hot core velocity, there 
was no obvious feature at the 2 \ - 1 1 frequency. Hence, CH3D was clearly not 

detected in this source. The upper limit to the column density and fractional 
abundance of CH3D on the Orion hot core, relative to H2, is N tot ~ 6 x 10 18 cm- 2 
and f ~ 6 x 10*6. These values imply that the methane abundance in the hot 
core is f < 6 x 10' 4 , consistent with current chemical models, which predict 
CH 4 /H 2 ~ 10-4. 

This search for CH3D also illustrates the need for careful and objective 

measurements in the detection of new interstellar molecules. Based on two lines 
alone (the Jk = lo - Oo and 2 q - lo transitions), we could have claimed an 
interesting discovery of interstellar CH3D. However, failure to observe a 
consistent third transition, namely the Jk = 2i - 1 1 line, showed beyond any doubt 
that such a claim would have been premature and absolutely incorrect. 

The authors also acknowledge grant NAGW -3065 from the NASA Origins 
of Solar Systems program. This research was also supported by funds from the 
Margaret Cullinan Wray Charitable Lead Annuity Trust. 
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Table 1: Lines Observed in the Bandpass 


Frequency (MHz) 

Identification 

Ta*(K) 

VlSR (km/ s ) AVi/ 2 (km/s) 

Comments 

229,647.7 

VyCN 

0.12 

5.3 

6.5 

LSB line; appears twice in 
spectrum 

U232,638 

— 

~ 0.1 

9.0 

~5 

blended with VyCN 

U232.647 

— 

0.1 

9.0 

3.8 

near CH 3 D; Jk= lo'Oo 
transition 

U465.232 

— 

0.5 

9.0 

~6 

blended with U 465,238 

U465.238 

— 

0.6 

9.0 

~6 

blended with U465,232 

U465,258 


0.5 

9.0 

5.8 

near CH 3 D; Jk= 2o-1q 
transition 

U468,120 



1.8 

9.0 

4.5 

— 


14 


Table 2: Observed Frequencies and Revised Constants for CH 3 D 


Transition (Jk) Frequency (MHz) 



lo-Oo 
2o- lo 

2 i - 1 1 


232,644.301 (0.075) 
465,250.691 (0.075) 
465,235.540 (0.075) 

Constant 


v (MHz)a) 

Previous Work (MHz) a ) 

Bo 

116,325.309(12) 

116,323(2)b), 116,325.308(27)c) 

Dj 


1.5725(21) 

1.572(7)b) 

Djk 


3.79(3) 

3.78(3)b) 


a) Errors listed are 3a. 

b) From Tarrago et al. (1976). 

c) From Pickett, Cohen, and Phillips (1980). 
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Figure Captions: 


Figure 1: Spectra taken at the frequencies of the CH 3 D transitions at 323 and 

465 GHz, respectively, using the CSO. The top panel shows the spectrum at 232 
GHz, where a spectral line occurs near the CH 3 DJK = lo - Oo frequency. A vinyl 

cyanide line is present in the LSB in these data, which appears twice due to an LO- 
shift. For one LO-setting, it is blended with a U-line from the USB. In the bottom 
panel, the data taken at 465 GHz are presented. Although a spectral feature occurs 
near the Jk= - lo frequency, there is not a line corresponding to the Jk = 2 i - 1 1 
transition. The frequency scale assumes Vlsr = 9 km/s in both cases. 

Figure 2: Laboratory spectrum of the CH 3 D:Jk= 2q - lo and Jk = ~ 1 1 

transitions near 465 GHz. These measurements confirm the 15 MHz splitting of 
the two K components of this transition. This spectrum represents one, 3 minute 
scan. 
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